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Abstract 


Marburg virus is the most deadly virus identified in the year 1967. MARV comprises a 19-KB 
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noninfectious single-stranded RNA genome that encodes seven basic structural proteins. Right 
now; there is no accessible treatment or vaccine to cure the MARV disease since couples of 
immunizations are under clinical trial studies. Here the computational approach is used to 
predict multi-epitope vaccine candidates against this virus. Both T cell and B cell were checked 
for the peptides to confirm that they can bring both humoral and cell-mediated immunity. A 


Omer T cell epitope ALSLTCAGI interacted with most of the MHC-I alleles and 10mer B cell 
epitope SLFVQAALYV was identified on antigenicity prediction. The regions; 7 — 12, 139 — 
144 and 186 — 191 amino acid residues were more accessible as B cell epitopes. The predicted 
epitope docked with the human HLA:*A0201. On the basis of the computational analysis, the 
predicted epitope can be the best accessible solution to cure infections of Marburg virus. 


Introduction 


The family Filoviridae incorporates two known genera, 
Ebola virus, and Marburg virus. The genus Ebola virus 
incorporates five species, while the Marburg virus genus 
incorporates single species: Marburg virus, which has two 
individuals: Marburg virus (MARV) and Ravn virus 
(RAVV) [1]. About Forty years ago, in August 1967, the 
first filovirus ever identified as the Marburg virus 
appeared in Europe, bringing about extreme and deadly 
hemorrhagic fever among the laborers of laboratory 
settings in Marburg and Frankfurt [2, 3]. Marburg 
hemorrhagic fever (MHF) is described by systemic viral 


replication, strange inflammatory reactions, and 
immunosuppression, including various systemic 
dysfunctions; hemorrhages, multi-organ failure, 


coagulation variations, edema, and shock, that often 
results in death [4]. 

MARV comprises a 19-KB noninfectious single-stranded 
RNA genome that encodes seven basic structural 
proteins: Nucleoprotein (NP), Virion protein 24 
(VP24),VP35, VP30, VP40, Glycoprotein (GP) and 
Polymerase L protein (L) [5]. Right now, there is no 
accessible treatment or vaccine to cure the MARV 
disease, whereas the couples of immunizations are under 


clinical trial studies. A novel method coordinating 
immunogenetics and immunogenomics with 
bioinformatics for the advancement of immunization 
vaccines is known as Vaccinomics [6]. This method is 
utilized for the discovery and development of new 
vaccines. The fast in-silico informatics-based method has 
received much recognition in the sequencing of numerous 
pathogenic genomes and protein sequence databases [7]. 
The "Vaccinomics" approach has sparingly turned out to 
be crucial for fighting infections, for example, multiple 
sclerosis [8], malaria [9], and tumors [10]. However, 
these techniques to develop vaccines work to determine 
the human leukocyte antigens (HLA) ligands and T cell 
epitopes [11] which determine the choice of the strong 
vaccine candidates concerned with the transporter of 
antigen presentation (TAP) molecules [12-15]. In the 
current study, an attempt is made to design an epitope- 
based peptide vaccine against the MARV strain of 
Marburg virus by means of the vaccinomics approach. 


Experimental 
Oany et al, [16] used novel and detailed method to design 


an epitope based vaccine against corona virus, in this 
research work similar methodology is applied to design an 
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epitope based vaccine against Marburg virus. The proper 
flowchart of applied method is shown in Figure 1. 


Selection of viral strain 

Viral Zone, a database of the ExPASy Bioinformatics 
Resource Portal was utilized to collect of MARV and its 
related data as well as the host, transmission, ailment, genus, 
family, genome, and proteome. The primary sequences of 
MARV proteins were retrieved from the UniProtKB 
database [17]. 


Determination of most antigenic protein 

VaxiJen v2.0 [18] a server for the forecast of defensive 
antigens and subunit vaccines, was utilized to identify the 
most powerful antigenic protein. 


T cell epitope Prediction 
The NetCTL 1.2 servers were utilized to determine the T 
cell epitopes [19]. 


Prediction of MHC1 binding alleles 

A tool of the Immune Epitope Database [20] was utilized to 
expect the MHC-I bindings. The stabilized matrix base 
method (SMM) [21] was utilized to find out the half- 
maximal inhibitory concentration (IC50) estimations of 
peptide bindings to MHC-I. For the binding examination, 
every one of the alleles was chosen, and the length was set at 
9.0. For the chosen epitopes, an online server was utilized to 
anticipate proteasomal cleavage, TAP transport, and MHC-I 
[22]. 
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Figure 1. The method applied to design a vaccine against 
Marburg Virus 


Allergen protein identification 

The AllerHunter server [23] was utilized to predict the 
allergenicity of proposed epitope to develop the vaccine. 
The Aller Hunter predicts allergens and additionally 
nonallergens with high specificity. This makes Aller Hunter 
is an exceptionally valuable tool for cross-reactivity 
prediction of allergens [24, 25]. 


B cell epitopes prediction 

The aim to predict the B cell epitope was to locate the 
potential antigen that would associate with B lymphocytes to 
start an immune response. Tools from IEDB were utilized to 
recognize the B cell antigenicity, including the Kolaskar and 
Tongaonkar antigenicity scale [26], Karplus and Schulz 
flexibility prediction [27], Emini surface accessibility 
prediction [28], The Chou and Fasman prediction tool [29] 
and Bepipred linear epitope prediction analysis [30]. 


Docking analysis of B cell and T cell epitopes 

For the docking studies, the T cell epitope ALSLTCAGI and 
B cell epitope SLFVQAALYV were subjected to PEP- 
FOLD electronic server [31] for 3D structural 
transformation, To investigate the interactions with various 
HLAs, a docking study was performed utilizing Galaxy 
PEPP-DOCK server. The HLA-A0201 was chosen for 
docking on the premise of the accessible Protein Data Bank 
(PDB) structure stored in the database, which interacted with 
proposed epitope [32]. 


Results and Discussion 


Results 

A total of seven primary sequences of NP, VP24,VP35, 
VP30, VP40, GP and L proteins were retrieved from 
UniProtKb database. The VaxiJen server evaluated all of the 
retrieved protein sequences to discover the most powerful 
antigenic protein. UniProtKB id:VP30 MARVM was 
selected as the most effective antigenic protein, with a 
maximum total prediction score of 0.5694. Then, this 
protein was utilized for further analysis. The NetCTL server 
anticipated the six effective T cell epitopes from the selected 
protein sequence such as NLGHILSYL, KLDETSLRA, 
ALSLTCAGI, SLINTMTEL, FVQAALYVM and 
VQASYDHFI with a specificity score of 0.970 and a 
sensitivity score of 0.80. The MHC-I alleles for which the 
epitopes showed higher affinity (IC50 <200 nM) are shown 
in Table 1. 

Amongst the six T cell epitopes, a 9mer epitope, 
ALSLTCAGI, interacted with most of the MHC-I alleles, 
including HLA:*H2Kb, HLA:*H2Kk, HLA:*A0203, 
HLA:* A0206, HLA:*A1101, HLA:* DRB0101, HLA:* 
DRB0701, HLA:*IAb, HLA:*IAs and TAP with higher 
affinity values of 0.75—1.0. The sequence-based allergenicity 
prediction tool specifically calculated by the use of the Aller 
Hunter tool, and the predicted queried epitope, allergenicity 
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score was 0.06 having a sensitivity of 98.4% and a 
specificity value of 13.1%. Different methods were used to 
predict potential B cell epitopes, the Kolaskar and 
Tongaonkar antigenicity prediction method examined 
antigenicity of B-cell epitope on the basis of the 
physiochemical properties of amino acids. The average 
antigenicity value obtained was 1.013. Nine epitopes 
represented to be a potential to express the B cell response 
shown in Table 2 and Figure 2. 


A 10mer B cell epitope SLFVQAALYV was designed on 
the basis of common amino acid residues produced in 
peptide sequences of Kolaskar and Tongaonkar antigenicity 
prediction. 

For a B cell to be a potential epitope, it must have proper 
surface accessibility. Therefore, Emini surface accessibility 
prediction was done. The regions from 7—12, 139-144 and 
186 — 191 amino acid residues were more accessible. The 
accessible residues are shown in Figure 3. 


Table 1. The six potential T-cell epitopes, along with their interacting HLA alleles and inhibitory concentration values 


Epitope Predicted IC50 Value (nM) Interacting HLA alleles 
NLGHILSYL 27.16, 53.06, 108.64, 179.06, 87.10, A0202, A0206, A0203, A0301 
10.79, 137.09 A1101, DRBO101, Iak 
KLDETSLRA 140.28, 83.18, 151.71, 1.58, 145.55, H2Kk, A1101, A6801, DRB0101, 
41.02, 88.72 DRB0701, IAd, IAs, 
ALSLTCAGI 39.99, 62.37, 9.73, 98.40, 62.23, H2Kb, H2Kk, A0203,A0206, A1101, 
46.13, 152.76, 186.64, 14.16, 8.30 DRBO101, DRB0701, IAb, IAs, TAP 
SLINTMTEL 217.77, 27.86, 48.98, 22.80, 128.82, A0202, A0203, A1101, DRB0101, 
110.92 IAs, TAP 
FVQAALYVM 43.55, 94.84, 21.73, 150.31, 1.16, A1101, A6801, A6802, B3501, 
45.39, 71.61, 34.99 DRBO101, DRB0401, IAd, IAs, 
VQASYDHFI 59.16, 146.89, 97.27, 83.37 A0202, A0203, A1101, Iak 


Table 2. B-cell epitopes based on Kolaskar and Tongaonkar antigenicity analysis 


Position Residue Start End Peptide 
249 A 246 252 VQAALYV 
247 Q 244 250 LFVQAAL 
159 S 156 162 HILSYLH 
248 A 245 251 FVQAALY 
245 F 242 248 ISLFVQA 
179 S 176 182 AALSLTC 
180 L 177 183 ALSLTCA 
244 L 241 247 SISLFVQ 
246 V 243 249 SLFVQAA 
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Figure 2. Kolashkar and Tongaonkar 
prediction of the most antigenic protein VP30 


antigenicity 


Chou and Fasman beta-turn prediction method predicted 
regions from 68 — 76 and 124 — 133 as potential Beta-turn 
regions. From the experimental confirmation, it is evaluated 
that the flexibility of the peptide is associated to its 
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Figure 3. Emini surface accessibility prediction of the most 


antigenic protein VP30 


antigenicity therefore, the Karplus and Schulz flexibility 
prediction method were used [39], this prediction method 
displayed the region of 273 — 280 amino acid residues as a 
most flexible region. 
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Table 3. Bepipred linear epitope prediction for B cell epitope 


St End Peptide 
1 16 MQQPRGRSRTRNHQVT 
18 20 TIY 
23 78 TQLPSKPHYTNYHPRARSMSSTRSSAESSPTNHIPRARPPSTFNLSKPPPPPKDMC 
87 99 CADPTCNRDHDLD 
115 146 LPNTDKTFRSPQDCGSPSLSKGLSKDKQEQTK 
167 169 GKL 
171 171 E 
205 214 NLPQDQNGVI 
224 224 D 
226 229 GGQF 
239 239 K 
257 257 P 
259 259 E 
261 263 SIS 
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Figure 4. Bepipred linear epitope prediction of the most 
antigenic protein VP30 


F 


Figure 6. Docking results of HLA:*A0201 and proposed B 
cell epitope SLFVQAALVY 


Bipred prediction method predicted that the peptide 
sequences from 68 — 76 and 124 — 133 amino acids are able 
to induce the preferred immune responses as B cell epitopes. 
These regions are shown in Table 3 and Figure 4. 


Figure 5. Docking results of HLA:*A0201 and proposed T 
cell epitope ALSLTCAGI. 


The T cell epitope ALSLTCAGI was selected on the basis 
of its interactions with large number of alleles, while the B 
cell epitope SLFVQAALYV was selected on the basis of 
frequent occurance of common amino acid residues, 
predicted as B cell epitope. The T cell epitope ALSLTCAGI 
and B cell epitope SLFVQAALYV docked with the human 
HLA:*A0201, they fixed into the alpha helix groove of 
allele properly that demonstrated their stability and efficacy 
as a potential vaccine against Marburg virus. The docked 
complexes are shown in Figure 5 and Figure 6. 

On the basis of all the computational analysis, it is suggested 
that the predicted epitope can be a best accessible solution to 
cure infections of Marburg virus. 


Discussion 

The development of novel vaccines in a shorter period of 
time is very important to cope up with the ever rising 
diseases [33-37]. Due to the advancement in the 
sequence-based technologies, now there is a lot of 
information available about the genomes and proteomes 
of several viruses. As a result, with the assistance of 
various bioinformatics tools, peptide-based vaccines can 
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be designed. Though, the concept of epitope-based 
vaccine design is getting recognition, but still, there is no 
work on epitope-based vaccine design for Marburg virus. 
In this work, an attempt is made to design Insilico 
epitope-based vaccine for Marburg virus. At present, 
vaccines are generally based B cell immunity. However, 
immunization taking into account T cell epitope has been 
energized, as the host can create a strong immune 
response by CD8+ T cell against the infected cell [38]. 
With time, because of antigenic drift, any foreign 
molecule can get away from the memory response of an 
antibody; however, the immune responses generated by T 
cell regularly give long lasting immunity. 

There are a few criteria that should be satisfied by an 
epitope-based vaccine candidate, and the proposed 
epitope against Marburg virus satisfied each of the 
criteria. Basically; allergenicity is one of the most 
important obstacles to the development of vaccines. 
Therefore, today, most of the vaccines designed to 
stimulate the immune system are checked for unfavorable 
allergic reaction, first [39] through inducing the T helper 
cell type 2 (Th2) and immunoglobulin E (IgE). 

The AllerHunter score is the likelihood that a specific 
sequence is a cross-receptive allergen. However, the limit 
for the prediction of allergen cross-reactivity is 
approximately greater than 0.08. Here, the proposed 
epitope's allergenicity score was 0.06, due to this, it was 
considered as a nonallergen. All these anticipated in 
silico results depend on the analysis of sequences and 
immune databases. This kind of study has gotten 
experimental validation [40] and hence, it is 
recommended that the proposed epitope would have the 
capacity to trigger a strong immune response as a peptide 
immunization in vivo. 


Conclusion 


In this research Insilico method was used to design the T cell 
and B cell epitope based peptide vaccine against Marburg 
virus. This study has revealed that combined computational 
approaches could be useful for predicting vaccine candidates 
against pathogens such as Marburg virus. In this way, in 
silico studies save both time and costs for development of 
vaccine candidates. They are helpful for researchers, and can 
guide the experimental work with higher possibilities of 
finding the desiredsolutions and with fewer clinical trials. In 
future this research work can be further utilized by the help 
of in-vitro experimental procetures to check its adequacy 
and efficacy 


Acknowledgement 

The Authors are grateful to Govt. Post Graduate College 
Mandian Abbottabad, Department of Bioinformatics and 
Bioinformatics Research Club of Government Post Graduate 
College Mandian Abbottabad for providing a platform to 
conduct research. 


Conflict of Interest 

This research work is unique and has not been submitted to 
any journal yet. None of the authors have challenged 
conflicts of interest. 


References 


1. Friedrich, B. M., Trefry, J. C., Biggins, J. E., Hensley, L. E., Honko, A. 
N., Smith, D. R., & Olinger, G. G.: Potential Vaccines and Post-Exposure 
Treatments for Filovirus Infections. Viruses 2012; 4(12): 1619-1650. 
http://doi.org/10.3390/v4091619 

2. Slenczka, W., & Klenk, H. D.: Forty Years of Marburg Virus. The Journal 
of Infectious Diseases 2007; 196(s2): S131-S135. 
http://doi.org/10.1086/520551 

3. Brauburger, K., Hume, A. J., Mühlberger, E., & Olejnik, J.: Forty-Five 
Years of Marburg Virus Research. Viruses 2012; 4(12): 1878-1927. 
http://doi.org/10.3390/v4101878 

4. Mehedi, M., Groseth, A., Feldmann, H., & Ebihara, H.: Clinical aspects 
of Marburg hemorrhagic fever. Future Virology 2011; 6(9): 1091-1106. 
http://doi.org/10.2217/fv1.11.79 

5. Thi, E. P., Mire, C. E., Ursic-Bedoya, R., Geisbert, J. B., Lee, A. C., 
Agans, K. N., et al.: Marburg virus infection in nonhuman primates: 
Therapeutic treatment by lipid-encapsulated siRNA. Science Translational 
Medicine 2014; 6(250):250ra116. 

6. Poland GA, Ovsyannikova IG, Jacobson RM.: Application of 
pharmacogenomics to vaccines. Pharmacogenomics 2009; 10(5):837—852. 

7. Flower DR: Bioinformatics for Vaccinology. Chichester: John Wiley & 
Sons, Ltd. 2008. 

8. Bourdette DN, Edmonds E, Smith C, et al.: A highly immunogenic 
trivalent T cell receptor peptide vaccine for multiple sclerosis. Mult Scler. 
2005; 11(5):552-561. 

9. Lopez JA, Weilenman C, Audran R, et al.: A synthetic malaria vaccine 
elicits a potent CD8(+) and CD4(+) T lymphocyte immune response in 
humans. Implications for vaccination strategies. Eur J Immunol. 2001; 
31(7):1989-1998. 

0. Knutson KL, Schiffman K, Disis ML.: Immunization with a HER-2/neu 
helper peptide vaccine generates HER-2/neu CD8 T-cell immunity in 
cancer patients. J Clin Invest. 2001; 107(4):477-484. 

1. Petrovsky N, Brusic V.: Computational immunology: The coming of age. 
Immunol Cell Biol. 2002; 80(3):248-254. 

2. Brusic V, Bajic VB, Petrovsky N.: Computational methods for prediction 
of T-cell epitopes — a framework for modelling, testing, and applications. 
Methods 2004; 34(4):436-443. 

3. Peters B, Bulik S, Tampe R, Van Endert PM, Holzhiitter HG.: Identifying 
MHC class I epitopes by predicting the TAP transport efficiency of 
epitope precursors. J Immunol. 2003; 171(4):1741-1749. 

4. Bhasin M, Raghava GP.: Analysis and prediction of affinity of TAP 
binding peptides using cascade SVM. Protein Sci. 2004; 13(3):596-607. 

5. Nielsen M, Lundegaard C, Lund O, Kesmir C.: The role of the 
proteasome in generating cytotoxic T-cell epitopes: insights obtained 
from improved predictions of proteasomal cleavage. Immunogenetics 
2005; 57(1-2):33-41. 

6. Oany, A., Pervin, T., & Emran, A.: Design of an epitope-based peptide 
vaccine against spike protein of human coronavirus: an in silico approach. 
Drug Design, Development and Therapy 2015; 1139. 
http://doi.org/10.2147/DDDT.S67861 

7. Apweiler R, Bairoch A, Wu CH, et al.: UniProt: the universal protein 
knowledgebase. Nucleic Acids Res. 2004; 32(Suppl 1):D115—D119. 

8. Doytchinova IA, Flower DR. VaxiJen: A server for prediction of 
protective antigens, tumour antigens and subunit vaccines. BMC 
Bioinformatics 2007; 8:4. 

9. Larsen MV, Lundegaard C, Lamberth K, Buus S, Lund O, Nielsen M.: 
Large-scale validation of methods for cytotoxic T-lymphocyte epitope 
prediction. BMC Bioinformatics 2007; 8:424. 

20. Buus S, Lauemegller SL, Worning P, et al.: Sensitive quantitative 
predictions of peptide-MHC binding by a ‘Query by Committee’ artificial 
neural network approach. Tissue Antigens 2003; 62(5):378-384. 

21. Peters B, Sette A.: Generating quantitative models describing the 
sequence specificity of biological processes with the stabilized matrix 
method. BMC Bioinformatics 2005; 6:132. 

22. Tenzer S, Peters B, Bulik S, et al.: Modeling the MHC class I pathway by 
combining predictions of proteasomal cleavage, TAP transport and MHC 
class I binding. Cell Mol Life Sci. 2005; 62:1025—1037. 


65 


23. 


24. 


25. 


26. 


27. 


28. 


29. 


30. 


Anum Munir et al., JIPBS, Vol 3 (4), 61-66, 2016 


Muh HC, Tong JC, Tammi MT.: AllerHunter: a SVM-pairwise system 
for assessment of allergenicity and allergic cross-reactivity in proteins. 
PLoS One 2009; 4(6):e5861. 

Liao L, Noble WS.: Combining pairwise sequence similarity and support 
vector machines for detecting remote protein evolutionary and structural 
relationships. J Comput Biol. 2003; 10(6):857—868. 

Chen YPP, Wong L, editors.: A Better Gap Penalty for Pairwise-SVM: 
Proceedings of the 3rd Asia-Pacific Bioinformatics, Singapore, January 
2005; 17-21. 

Thévenet P, Shen Y, Maupetit J, Guyon F, Derreumaux P, Tufféry P.: 
PEP-FOLD: an updated de novo structure prediction server for both linear 
and disulfide bonded cyclic peptides. Nucleic Acids Res.2012; 40(Web 
Server issue): W288-—W293. 

Berman HM, Westbrook J, Feng Z, et al.: The protein data bank. Nucleic 
Acids Res. 2000; 28(1):235—242. 

Kolaskar AS, Tongaonkar PC.: A semi-empirical method for prediction of 
antigenic determinants on protein antigens. FEBS Lett. 1990; 276(1— 
2):172-174. 

Emini EA, Hughes JV, Perlow DS, Boger J.: Induction of hepatitis A 
virus-neutralizing antibody by a virus-specific synthetic peptide. J Virol. 
1985; 55(3):836-839. 

Karplus PA, Schulz GE.: Prediction of chain flexibility in proteins. 
Naturwissenschaften. 1985; 72(4):212-213. 


31. 


39. 


40. 


Larsen JE, Lund O, Nielsen M.: Improved method for predicting linear B- 
cell epitopes. Immunome Res. 2006; 2(1):2. 

Chou PY, Fasman GD.: Empirical predictions of protein conformation. 
Annu Rev Biochem. 1978; 47:251-276. 

Marshall SJ.: Developing countries face double burden of disease. Bull 
World Health Organ 2004; 82(7):556. 

De Groot AS: Rappuoli R.Genome-derived vaccines. Expert Rev 
Vaccines 2004; 3(1):59-76. 

Korber B, LaBute M, Yusim K: Immunoinformatics comes of age. PLoS 
Comput Biol. 2006; 2(6):e71. 

Fauci AS.: Emerging and re-emerging infectious diseases: influenza as a 
prototype of the host-pathogen balancing act. Cell 2006; 124(4):665-670. 
Purcell AW, McCluskey J, Rossjohn J.: More than one reason to rethink 
the use of peptides in vaccine design. Nat Rev Drug Discov. 2007; 6(5): 
404-414. 

Shrestha B, Diamond MS.: Role of CD8+ T cells in control of West Nile 
virus infection. J Virol. 2005; 78(15):83 12-8321. 

McKeever TM, Lewis SA, Smith C, Hubbard R.: Vaccination and allergic 
disease: a birth cohort study. Am J Public Health 2004; 94(6):985—989. 
Khan MK, Zaman S, Chakraborty S, et al: In silico predicted 
mycobacterial epitope elicits in vitro T-cell responses. Mol Immunol. 
2014; 61(1):16—22. 


66 


